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ABSTRACT 

In  this  report  we  develop  alternative  wholesale  provisioning  models  for 
secondary  items  with  an  aim  toward  use  by  both  the  Ships  Parts  Control  Center 
(SPCC)  and  the  Aviation  Supply  Office  (ASO).   General  background  into  the 
wholesale  provisioning  problem  is  given,  and  the  present  provisioning  models 
in  use  at  SPCC  are  summarized.   We  review  the  demand  forecasting  problem,  and 
we  address  the  issue  of  demand  distribution  assumptions.   Then  various 
alternative  wholesale  provisioning  models  are  developed.    These  models 
attempt  to  allocate  optimally  a  given  provisioning  budget  with  such  objective 
functions  as  (1)  units  short,  (2)  essentiality-weighted  units  short,  (3) 
essentiality-weighted  requisitions  short,  (4)  supply  material  availability 
(SMA),  (5)  time-weighted  units  short,  (6)  essentiality-weighted  mean  supply 
response  time,  and  (7)  system  availability.   We  conclude  with  theoretical 
comments  about  the  alternative  models  and  discussion  of  a  companion  problem 
that  can  be  solved  to  determine  the  minimum  budget  required  to  attain  a 
specified  level  of  performance. 

This  is  the  first  of  three  reports  on  the  wholesale  provisioning  problem. 
The  second  report  will  present  results  of  a  numerical  evaluation  of  the 
proposed  alternative  models  and  the  current  model.   The  thira  report  addresses 
the  issue  of  phased  provisioning  and  offers  a  heuristic  method  for  determining 
how  to  phase  budget  expenditures  over  time. 


1 .   INTRODUCTION 

In  the  spring  of  1982,  the  Naval  Supply  Systems  Command  (NAVSUP)  asked 
the  Naval  Postgraduate  School  to  develop  improvements  to  the  existing 
peacetime  wholesale  provisioning  models  for  secondary  items  used  by  the  Ships 
Parts  Control  Center  (SPCC)  and  the  Aviation  Supply  Office  (ASO) .   This  effort  was 
motivated  by  NAVSUP' s  Resolicitation  Project,  the  major  objective  of  which 
is  to  acquire  new  computer  hardware.   However,  it  also  provides  an  opportunity 
to  take  a  hard  look  at  existing  models  in  use  by  NAVSUP 's  Inventory  Control 
Points  (ICP)  at  SPCC  and  ASO  and  to  make  changes  as  desired  before  the  ICP 
software  is  installed  in  the  new  computers.   These  changes  should  include 
wholesale  and  retail  model  improvements,  both  in  provisioning  and 
replenishment  of  secondary  items  .   The  most  important  questions  to  be 
answered  by  the  Naval  Postgraduate  School  are: 

1.  Do  the  existing  models  provide  the  best  supply  support  for  the  .lavy '  s 
weapon  systems? 

2.  How  should  those  models  which  do  not  provide  the  best  support  be 
changed? 

In  responding  to  these  questions,  it  makes  sense  to  begin  with  "the 
beginning"  of  supply  support —  provisioning  for  a  weapon  system — because,  once 
an  improved  provisioning  model  is  developed,  the  replenishment  models  can  be 
modified  in  such  a  way  as  to  provide  a  smooth  transition  to  continued 
wholesale  supply  support.   This  report  will  therefore  concentrate  on  the 
provisioning  models. 

Wholesale  supply  support  provides  backup  or  "systems"  stock  of  repairable 
and  consumable  items  for  the  retail  levels.   At  present  the  retail  levels  use 
models  such  as  CGSALs  (Coordinated  Shipboard  Allowance  List)  and  A7CALS 
(Aviation  Coordinated  Allowance  List).   These  are  stocks  of  3pare  secondary 


items  (items  which  may  tail  in  use)  designed  to  provide  sufficient  support  for 
a  ship  deployment  period  of,  say,  90  days.   Whenever  an  item  is  used  by  a  ship 
or  a  squadron,  a  requisition  is  immediately  submitted  to  the  wholesale  system 
for  a  replacement.   Thus,  the  COSALs  and  AVCALs  serve  only  as  emergency 
protection.   All  demands  ultimately  must  be  satisfied  by  the  wholesale  stocks. 

The  wholesale  provisioning  model  provides  protection  until  the  wholesale 
replenishment  model  can  initiate  the  first  replenishment  buy  (or  repair 
quantity  for  those  secondary  items  which  can  be  repaired)  and  that  buy  is 
manufactured  and  delivered.   Thus,  wholesale  provisioning  buy  protection 
begins  on  the  date  of  preliminary  operational  capability  (POC)  and  continues 
for  at  least  a  replenishment  procurement  (or  repair)  leadtime. 

Actually,  it  is  expected  that  a  wholesale  replenishment  buy  will  be 
initiated  some  substantial  time  after  POC.   This  will  occur  when  the  reorder 
point  established  by  the  replenishment  model  is  reached.   The  initial  reorder 
point  is  determined  differently  by  the  two  ICP's  and  it  may  be  based  on  the 
initial  expected  demand  which  is  predicted  as  part  of  the  provisioning 
process.   However,  after  six  months,  the  reorder  point  begins  to  reflect  the 
actual  demands  placed  on  the  supply  system.   This  may  result  in  the  value  of 
the  reorder  point  decreasing  for  a  year  or  two. 

The  provisioning  buy  must  be  made  early  enough  so  that  the  stock  will  be 
on  hand  by  the  POC  date  or  shortly  thereafter.   Thus  it  must  be  made  at  least 
a  procurement  leadtime  prior  to  that  date.   Typical  values  of  this  initial 
leadtime  appear  to  be  between  2.0  and  2.5  years.   With  such  long  leadtimes  the 
quality  of  forecasted  installation  schedules  of  a  new  system  may  be  poor. 
And,  when  protection  must  be  provided  for  a  subsequent  time  which  is  at  least 
as  long  as  a  replenishment  leadtime,  the  quality  of  the  forecasted  demands 
over  that  time  can  also  be  ooor. 


Provisioning  guidance  is  provided  by  the  Department  of  Defense 
Instruction  (DODINST)  4140.42  [7].   The  authors  of  [7]  were  of  the  opinion 
that  much  stock  of  secondary  items  was  procured  which  was  not  subsequently 
needed.   Thus,  in  an  attempt  to  remedy  this  perceived  situation,  a  model  was 
developed  which  would  provide  only  a  very  conservative  quantity  of  wholesale 
stock  for  those  items  which  were  new  to  the  supply  system.   For  items  already 
stocked  in  the  system,  they  allowed  no  additional  stock  to  be  procured 
specifically  for  a  new  system.   However,  if  the  anticipated  increase  in  demand 
was  sufficiently  large,  it  could  be  entered  into  the  demand  forecasting 
process  by  the  inventory  manager  and  its  influence  could  trigger  an  earlier 
replenishment  buy  than  would  have  occurred  without  it. 

The  authors  of  [7]  did  appreciate  that  improved  models  for  provisioning 
new  items  would  evolve  and  stated  that: 

Inventory  optimization  models  may  be  used  as  a  basis  for  stockage  in  lieu 
of  the  model  in  this  Instruction  if  the  following  criteria  are  met: 
(1)  An  optimization  technique,  developed  to  minimize  system  downtime  or 
time-weighted  requisitions  short  is  used.   (2)  No  lower  limits  are  placed 
on  requirements  that  will  result  in  stockage  as  a  aemand-based  item,  of 
an  item  that  would  not  be  stocked  without  the  lower  limits.   (3)  A 
monetary  base  point  is  determined  indicating  the  approximate  value  of  the 
requirement  in  accordance  with  the  instructions  contained  in  enclosure  2. 
The  second  criterion  may  be  omitted  if  only  those  items  meeting  the 
stockage  criteria  in  accordance  with  this  Instruction  are  considered  in 
the  optimization  model.  The  longest  applicable  Program  F orecast  Period 
may  be  used  in  the  optimization  model,  after  the  monetary  base  has  been 
es tablishea. 


Such  improvements  have  indeed  been  developed  by  NAVSUP  and  are  in  use  at 
both  SPCC  and  ASO.   However,  as  was  stated  above,  NAVSUP  is  interested  in 
developing  even  better  models.   This  report  represents  the  first  step  and  will 
concentrate  on  a  model  to  improve  the  initial  supply  support  for  new  secondary 
items  that  will  be  needed.   The  issue  of  how  to  handle  items  already  stocked 
will  be  the  subject  of  a  later  report. 

This  report  begins  with  a  discussion  of  the  existing  wholesale 
provisioning  model  used  by  SPCC  since  it  is  more  mathematically  sophisticated 
than  the  ASO  model  and  it  is  also  easier  to  understand.   Next,  the  demand 
process  is  examined  and  several  alternative  models  intended  for  use  by  both 
SPCC  and  ASO  are  developed.   These  alternative  models  develop  a  budget  using 
the  COSDIF  process  and  then  allocate  that  budget  so  as  to  optimize  system 
performance  with  respect  to  such  measures  as  supply  material  availability, 
mean  supply  response  time,  and  availability.   A  discussion  of  some  of  the 
theoretical  properties  of  the  models  is  provided  and  recommendations  based  on 
these  properties  are  made. 

Three  appendices  are  included.   The  first  presents  the  details  of  the 
forecasting  models  from  [7]  and  from  the  SPCC  and  ASO  models.   The  second 
discusses  the  cost  difference  or  COSDIF  formula  which  is  required  by  [7]  for 
developing  the  budget  constraint.   The  last  examines  a  widely  held 
misconception  about  the  relationship  between  mean  supply  response  time  and 
availability. 

This  is  the  first  of  three  reports  en  wholesale  provisioning.   The  second 
report  describes  the  computational  algorithms  used  for  provisioning  with  each 
of  the  alternative  models  and  provides  some  numerical  evaluations  of  the 
models.   The  third  report  examines  the  issue  of  phased  provisioning  and  offers 
a  heuristic  procedure  for  determining  how  to  divide  a  provisioning  budget  over 
a  period  of  time. 


2.   THE  CURRENT  SPGC  PROVISIONING  MODEL 
DEVELOPMENT  OF  THE  BUDGET  CONSTRAINT 

The  model  that  SPCC  uses  to  determine  initial  wholesale  stocks  for  new 
items  is  a  variant  of  the  model  proposed  in  reference  [7],   The  application/ 
operation  which  implements  the  model  is  designated  as  D55  [19].   The  model 
makes  range  and  depth  decisions  separately  for  each  item  in  the  system  for  the 
purpose  of  developing  a  budget  constraint. 

The  initial  phase  of  the  process  is  to  develop  a  budget  constraint.   The 
first  step  of  this  phase  is  the  determination  of  a  range  of  items  to  be 
considered  for  stockage.   This  is  done  by  computing  the  COSDIF  value  for  each 
item.   The  purpose  of  the  computation  is  to  decide  if  an  item  is  to  be  stocked 
as  a  demand-based  item.   The  COSDIF  value  is  obtained  from  a  formula  which  was 
stated  in  reference  [7],   The  current  version  of  COSDIF  used  by  SPCC  is 
presented  in  reference  [19]  and  is  reproduced  in  Appendix  B  of  this  report. 
Appendix  B  also  presents  the  motivation  for  the  formula  and  a  critique  of  the 
procedure.   Elements  of  the  COSDIF  include  the  expected  costs  of  stocking  the 
item  and  having  no  demands  over  two  years,  the  expected  costs  of  stocking  an 
item  and  experiencing  demand  for  one  year,  and  the  expected  costs  of  not 
stocking  an  item  for  one  year  and  having  to  make  spot  buys  when  demands  occur. 
The  difference  between  the  sum  of  the  first  two  elements  and  the  last  element 
is  the  "cost  difference"  or  COSDIF.   If  the  value  of  the  COSDIF  is  negative 
then  the  item  will  be  stocked;  otherwise  it  will  not  be  stocked  insofar  as 
budget  development  is  concerned. 

If  an  item  has  been  determined  by  the  COSDIF  as  appropriate  to  stock, 
its  depth  of  stockage  is  next  computed.   The  depth  is  taken  to  be  equal  to  the 


demand  expected  over  a  period  of  time  consisting  of  the  forecasted  replenish- 
ment procurement  ieadtime  plus  one  quarter.   The  extra  quarter  of  demand  is 
viewed  as  safety  stock.   The  procedure  used  to  develop  the  initial  demand 
prediction  is  discussed  in  Appendix  A. 

Those  items  which  fail  to  satisfy  the  C0SD1F  criterion  are  next 
re-examined  to  determine  if  they  are  coded  as  insurance  or  numeric 
stockage  objective  (NSO)  items.   If  so,  the  depth  of  each  item  is  taken  to  be 
one  minimum  replaceable  unit  (MRU). 

The  cost  of  buying  the  amount  of  each  item  to  the  depth  specified  above 
is  then  determined,  and  the  total  dollar  value  of  a  provisioning  package  is 
obtained  by  summing  these  procurement  costs  over  all  of  the  NSO  items,  the 
insurance  items,  and  those  items  passing  the  COSDLF  test.   This  sum  is  then 
established  as  "the  budget  constraint"  in  accordance  with  reference  [7],   The 
dollar  value  so  determined  is  considered  as  a  firm  upper  limit  on  the  amount 
of  money  that  can  be  spent  to  purchase  those  items  that  are  to  be  stocked. 

The  actual  range  and  depth  of  items  that  must  be  stocked,  however,  may 
deviate  fron  the  values  used  to  determine  the  budget.   As  was  noted  in  chapter 
I,  an  alternative  provisioning  model  may  be  used  to  implement  the  range  and 
depth  provided  the  budget  constraint  as  determined  above  is  satisfied  and 
the  model  objective  is  consistent  with  minimization  of  system  down  time  or 
time-weighted  requisitions  short.   The  determination  of  the  budget  constraint 
is  therefore  an  important  element. 

THE  VARIABLE  THRESHOLD  MODEL  FOR  BUDGET  IMPLEMENTATION 

SPCC  has  developed  a  model,  called  the  Variable  Threshold  Model,  for 
making  the  budget  implementation.   In  that  model,  all  new  items  are  considered 


to  be  candidates  for  stockage  and  are  re-evaluated  based  on  unit  cost  and  the 
expected  demand  during  the  procurement  leadtime. 

The  first  step  of  the  model  is  to  compute  the  variable  threshold  value  of 
each  item.   This  value  is  computed  using  the  variable  threshold  formula: 

V(i)  =   (l  -  exp(-Di))  /  C. 

where  C.  is  the  unit  cost  of  the  item  and  D-j_  is  the  expected  demand  during  the 
procurement  leadtime  (see  Appendix  A).   The  variable  threshold  value  is  the 
mechanism  used  by  SPCC  to  rank  the  items  in  a  provisioning  package  according 
to  desirability  of  stockage.   Those  items  having  a  high  value  are  placed  at 
the  top  of  the  list.  The  numerator  of  the  threshold  value  is  the  probability 
of  one  or  more  demands  during  a  procurement  leadtime  when  the  demand  is 
Poisson  distributed.   It  comes  from  an  analysis  by  Burton  and  Jaquette  [2]  in 
which  the  term  represented  the  decrease  in  the  expected  number  of  backorders 
at  a  random  time  resulting  from  an  increase  in  units  stocked  from  zero  to  one 
for  an  item  that  has  Poisson  distributed  demands.   Reference  [2]  addressed 
only  very  low-demand  items  which  are  currently  viewed  as  insurance  or  NSO 
items. 

The  variable  threshold  formula  is  a  typical  "benefit-to-cost"  ratio  used 
in  the  first  step  of  a  marginal  analysis  allocation  scheme.   Similar 
expressions  are  found  in  the  alternative  models  which  we  present  later. 

The  next  step  is  to  compute  the  depth  for  each  item  based  on  a  risk 
formula  similar  to  that  used  in  the  replenishment  models.   The  formula  is 


RISK  =  IC  /  (IC  +  XSE) 


where  IC  is  the  annual  holding  cost,  \s    is  the  shortage  cost,  and  E  is  th< 
item  essentialitv. 


Let  H(x)  represent  the  probability  that  the  leadtime  demand  for  a  given 
item  is  greater  than  or  equal  to  x.   Then,  the  depth  for  the  item  is  taken  to 
be  the  smallest  value  of  x  such  that  H(x)  is  not  greater  than  the  calculated 
RISK  value.   Thus,  if   x.*   represents  the  depth  for  item   i, 


x.  "  =  mm  x; 

1 


IC. 
such  that  H.(x)  <  RISK  = 


i    -        IC.  +  A.E. 
1    li 


The  actual  probability  distribution  used  for  the  calculation  depends  on 
the  expected  annual  demand  frequency.   The  Poisson  distribution  is  used  if  the 
expected  annual  demand  is  less  than  or  equal  to  one;  the  normal  distribution 
is  used  when  the  expected  annual  demand  or  procurement  leadtime  demand  is  20 
or  more;  and  the  negative  binomial  distribution  is  used  when  the  expected 
annual  demand  is  between  one  and  20.   The  actual  depth  is  constrained  to  be 
no  more  than  two  years  of  expected  demand  if  the  item  is  consumable  and  no 
more  than  the  expected  "net"  attrition  demand  during  a  procurement  leadtime 
plus  one  quarter  if  the  item  is  repairable  (see  Appendix  A) . 

The  final  step  is  to  allocate  the  budget  constraint  among  the  items.   The 
part  of  the  total  budget  required  to  purchase  the  MRU's  for  the  insurance  and 
NSO  items  is  reserved  for  these  items  and  is  subtracted  from  the  total  budget 
available.   The  remaining  dollars  are  next  allocated  by  successively 
purchasing  the  depths  determined  above  for  the  items  ranked  in  order  of 
highest  to  lowest  variable  threshold  value.   The  item  with  the  highest 
variable  threshold  value  is  purchased  in  the  depth  specified  above;  the 
remaining  budget  is  decremented;  the  item  with  the  second  highest  variable 
threshold  value  is  purchased  next.   This  procedure  continues  until  the  budget 
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is  depleted  or  until  every  item  is  purchased  at  the  specified  depth.   Ai 
items  already  designated  as  insurance  or  NSO  items  that  qualify  for  stockage 
under  the  above  procedure  will  first  have  the  depth  quantities  reduced  by  the 
already  funded  MRU's. 

A  comparison  of  the  SPCC  budget  allocation  approach  with  that  described 
in  reference  [7]  was  done  by  the  Fleet  Material  Support  Office  (FMSO)  and 
reported  in  reference  [9],   The  results  of  the  study  showed  that  the  variable 
threshold  model  provided  more  range  and  less  depth  than  did  the  model  of 
reference  [7]  for  the  same  budget  amount. 


3.   MODEL  BACKGROUND 

INTRODUCTION 

The  provisioning  model  currently  in  use  at  SPCC  was  discussed  in  the 
previous  chapter.   That  model  is  a  variant  of  the  model  proposed  in  DODI 
4140.42.   The  provisioning  procedure  described  in  DODI  4140.42  is  an  aggregate 
of  the  Army's  replenishment  range  model  (COSDIF)  and  an  attempt  to 
conservatively  forecast  leadtime  demand  (under  the  assumptions  that  the 
components  have  reliabilities  equal  to  their  design  values).   That  instruction 
provides  uniform  guidelines  for  the  military  services  under  which  the  services 
could  develop  their  own  provisioning  models  for  carrying  out  the  budget 
implementation  The  instruction  does,  however,  allow  only  certain  types  of 
models.   Paragraph  IV. D  of  the  instruction  states: 

"Mathematical  models  may  be  used  which  may  provide  for  a  different  mix  of 
inventory  than  that  stated  herein,  providing  a  financial  base  is 
established  in  accordance  with  the  policies  outlined  in  this  instruction 
and  an  objective  to  minimize  system  downtime  or  time-weighted 
requisitions  short  is  included  in  model." 

We  examine  models  in  this  report  which  focus  on  these  types  of  objectives.   In 
particular,  we  develop  models  which  optimize  with  respect  to: 

1 .  expected  number  of  stockouts 

2.  expected  time-weighted  units  (requisitions)  short 

3.  mean  supply  response  time,  and 

4.  availability. 

Before  presenting  these  models,  it  is  appropriate  to  examine  some  background 
issues  which  affect  all  of  the  models  that  follow. 
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THE  PROVISIONING  INTERVAL 

There  are  four  key  dates  for  consideration  in  the  provisioning  problem: 

1.  POC:   the  date  of  planned  operational  capability  (part  of  the 
provisioning  package  should  be  in  place  by  this  date), 

2.  POC  -  PCLT:   a  procurement  leadtime  prior  to  PCC , 

3.  POC  +  TR:   the  time  beyond  POC  at  which  the  first  replenishment  buy 
is  made ,  and 

4.  POC  +  TR  +  PCLT:   the  time  at  which  the  first  replenishment  buy  is 
received. 

The  primary  concern  of  provisioning  is  to  determine  the  number  of  units 
of  each  component  of  a  system  to  be  purchased  at  time  POC  -  PCLT  (and  possibly 
purchase  at  additional  times  between  POC  -  PCLT  and  POC).   In  order  to  reduce 
the  likelihood  of  stockouts,  the  provisioning  buy  made  at  time  POC  -  PCLT 
should  be  sufficiently  large  to  satisfy  anticipated  demands  in  the  interval 
(POC,  POC  +  TR  +  PCLT).   If  good  forecasts  of  demand  during  that  interval  and 
of  the  times  TR  and  PCLT  were  available,  the  solution  of  the  provisioning 
problem  would  be  fairly  routine.   However,  the  provisioning  problem  is 
characterized  by  a  great  deal  of  uncertainty  with  respect  to  the  failure  rates 
(BRF'S)  of  new  equipments,  the  procurement  leadtime,  the  time  TR ,  and  even  the 
installation  schedule  of  new  equipments. 

During  the  time  interval  from  POC  to  POC  ~   TR  +  PCLT,  the  population  of 
equipments  will  grow  as  additional  units  are  installed  in  the  operational 
community.   This  will  cause  the  aggregate  failure  rates  to  increase  over  the 
time  interval  of  concern  in  the  provisioning  problem.   Actual  installation 
schedules  are  subject  to  a  great  deal  of  uncertainty;  the  3RF  estimates  ofcen 
represent  only  guesses  at  the  failure  rates  since  there  are  little  or  no 
operational  data  from  which  to  estimate  the  failure  rates;  and  the  time  lags 
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are  very  large  (the  time  interval  from  POC  -  PCLT  to  POC  t  TR+  PLT  may  be  on 
the  order  of  four  years  or  more).   Thus  it  is  not  surprising  that  the 
provisioned  quantities  frequently  don't  match  the  demands  very  well.  Faced 
with  such  uncertainties,  no  analytical  model  will  cure  all  of  the  provisioning 
ills.   However,  the  use  of  rational  provisioning  procedures  should  minimize 
such  problems  and,  in  addition,  lead  to  improved  replenishment  models. 

PROVISIONING  OBJECTIVE  FUNCTIONS 

A  key  ingredient  of  DODI  4140.42  is  the  budget  constraint  for  the 
procurement  of  items  to  be  provisioned.   No  flexibility  exists  in  the  method 
used  to  determine  the  budget.   (The  actual  values  of  the  parameters  used  in 
the  COSDIF  formula  can,  however,  be  selected  by  the  users.)   As  was  noted  at 
the  start  of  this  chapter  DODI  4140.42  does  allow  limited  flexibility  in  the 
type  of  provisioning  model  used  to  allocate  the  given  budget  among  new  items 
being  considered.   It  can  either  minimize  time-weighted  units  short  or 
minimize  system  downtime.  However,  models  such  as  the  variable  threshold  which 
actually  do  neither  have  been  approved.   As  a  consequence,  many  measures  of 
system  performance  have  been  and  continue  to  be  considered  for  objectives  of 
provisioning.   Some  of  the  more  popular  are  availability,  probability  of  a 
stockout,  fill  rate,  backorder  accumulation,  system  downtime,  time-weighted 
requisitions  short,  essentiality-weighted  shortages,  mean  supply  response 
time,  and  inventory  costs.   The  use  of  costs  as  the  principle  measure  of 
performance  is  complicated  by  the  impossible  task  of  quantifying  such  elements 
as  stockout  costs  and  obsolescence  costs.   The  use  of  true  system  availability 
models  for  the  wholesale  level  is  complicated  be  the  need  for  information  at 
the  operational  level  on  a  system-by-system  basis.   That  information  is 
currently  not  available.   Such  models  are  therefore  not  presently  feasible  for 
the  wholesale  level. 


.2 


Some  implementations  of  availability  models  for  stockage  decisions  define 
availability  for  individual  components  (as  opposed  to  weapon  systems)  in  terms 
of  the  ratio  of  mean  time  between  failures  (MTBF )  to  the  sum  of  MTBF  ,  mean 
time  to  repair  (MTTR) ,  and  mean  supply  response  time  (MSRT)  (or  other  similar 
expressions);  see,  for  example,  references  [  3]  and  [23]  .   The  proponents  of 
such  availability  models  argue  that  MSRT  is  the  only  term  which  pertains 
directly  to  the  provisioning  stockage  decision.   Consequently,  these 
availability  models  really  attempt  to  minimize  MSRT.   (Appendix  C  provides  a 
simple  example  to  show  that  minimization  of  MSRT  does  not  necessarily  yield 
the  same  solution  as  maximization  of  system  availability.) 

FORM  OF  THE  MODELS 

In  all  of  the  models  which  we  develop,  the  optimization  problem  is  of 

the  form: 

minimize    (maximize)  V      f.(s.) 

.11 

l 

subject    to  T      c.s.      <   B. 

J  L.        i  i     — 

i 

where  fi(si)  is  the  performance  level  for  item  i  when  s\_   units  are  stocked,  c^ 
is  the  unit  price  of  item  i,  and  3  is  the  total  amount  budgeted  for  the 
provisioning  package.  The  form  of  the  optimization  problem  implicitly  assumes 
separability;  that  is,  the  total  system  performance  can  be  separated  into  a 
function  which  is  the  sum  of  the  performance  generated  by  each  of  its 
components.   This  form  may  not  always  be  appropriate  because  it  ignores  the 
system  configuration  and  any  interaction  between  components.   However,  we  feel 
that  such  an  assumption  is  necessary  to  assure  mathematical  tractibility  and 
implementation  feasibility.   It  is  also  in  keeping  with  the  piece-parts 


[2 


support  view  which  is  traditional  in  the  supply  system.   True  "systems 
support"  models  are  still  a  long  time  away. 

When  the  performance  functions  are  separable  and  well  behaved  the 
mathematical  solution  to  the  optimization  problem  can  readily  be  obtained 
using  such  straightforward  techniques  as  the  generalized  Lagrange  multiplier 
approach  of  Everett  [  8]  or  marginal  analysis  [12]. 

RANGE  RULES 

Ail  of  the  models  which  we  will  address  solve  simultaneously  the  range 
and  depth  problems.   There  is  no  explicit  range  rule.   Instead,  zero  depth  is 
equivalent  to  a  no-stockage  decision.   Insurance  items  and  numeric  stockage 
objective  (NSO)  items  will  have  to  be  handled  separately  from  demand  based 
items  just  as  they  are  today.  For  the  items  which  are  demand  based  the  models 
require  an  assumption  about  the  distribution  of  demand.    The  next  chapter 
supports  the  use  of  the  Poisson  probability  distribution  with  an  expected 
demand  rate  which  depends  on  the  BRF  estimate  and  a  time-weighted  installed 
population  (TWAMP). 
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4.   THE  DEMAND  DISTRIBUTION 

One  of  the  most  critical  elements  in  any  stockage  model  is  the  estimation 
of  the  total  number  of  failures  (demands)  for  a  component  which  will  be 
generated  during  the  provisioning  protection  interval.   The  proper  way  to 
handle  demand  uncertainty  is  to  model  the  demand  process  as  a  random  process 
with  a  probability  law  from  a  known  family  of  probability  distributions  and 
having  parameters  which  are  estimated  from  the  specific  information  that  is 
available  for  each  item.   In  the  provisioning  problem,  this  information  is 
pretty  skimpy,  consisting  essentially  of  contractor  provided  estimates  of 
failure  rates  (TRF's)  and  the  schedule  of  installations. 

This  chapter  discusses  the  probability  distribution  that  should  be  used 
for  the  provisioning  stockage  decisions.    We  show  that,  under  the  assumption 
of  a  constant  failure  rate,  the  total  demand  for  a  given  item  over  any  time 
interval  is  Poisson  distributed  with  parameter  A  which  depends  on  the 
time-weighted  average  month's  program  (TV/AMP). 

Consider  a  single  component  with  constant  failure  rate  a  during  the 
interval  (0,T].   Let  N  be  the  total  number  of  installations  of  the  component 
during  (0,T],  and  let  0  =  t\    <  t2  <  • ♦ •  <  tft  be  the  installation  times  of  the 
N  units  of  the  component.   Let  X^  be  the  random  number  of  failures  of  unit  i 
in  (0,T].   Then,  because  of  the  constant  failure  rate,  the  distribution  of  Kj 
is  Poisson  with  parameter  ct(T-t^)  (see,  for  example,  Ross  [21]);  i.e., 

P(Xi  =  x)  =  (a(T-ti))xexp(-a(T-ti))  /  x!        x=l,2,... 

Let  us  now  consider  the  total  number  of  demands  generated  by  the  i^  and 
j th  installed  units.  If  we  assume  that  the  failures  of  the  ita  and  2  units 
are  independent,  the  probability  distribution  of  the  sum  of  Xj_  and  Xj  is  given 
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by  the  convolution  of  the  distributions  of  X^  and  Xj .   Thus 


m 


P(X.   +  X.   =  m)  =   I   P(X.  =  m-k)P(X.  =  k) 
1     J  i,=n   J  i 


k=0 


m         m-k  k 

J(o(T-tj))    exp(-a(T-tj))  (a(T-ti))  exp(-a(T-ti) ) 
k=0- 


(m-k)! 


v  i 


.m-k 

m   ,  /  l— t 


am(2T-t.-t.)mexp(-»(2T-t.-t.))l(°)' 


V^i  V 


m!  k=0   \  2T-t.-t.l   \2T-t.-t. 

1  J   '     i   J 


If  we  let  p  =  (T-ti)  /  (2T-ti-tj),  then  1-p  =  (T-cj)  /  (2T-tj-tj)  and  the 

above  summation  is  recognized  to  be  the  sum  of  binomial  probability  masses. 

Hence,  that  sum  is  unity,  and  the  distribution  of  the  sum  of  X^  and  Xj  is 

Poisson  with  parameter  a(2T-t^-tj).   Let  D  =  X^  +  X2  +  . .  .  +  Xjj  be  the  total 

demand  over  all  installed  units  during  the  interval  (0,T],   The  above  result 

is  easily  extended  to  show  that  D  is  Poisson  distributed  with  parameter 

m 
a(NT-  ]  t^).   Let  f(t)  be  the  total  number  of  installations  at  time  t.   Then 
1-1 

f(t)  is  a  step  function  with  jumps  at  times  t^,  to,  •••,  tjsj .   The  integral  of 

N  N 

f(t)  over  the  interval  is  therefore    (T-t^)  =  NT  -  V  tj_.   But  this  integral 

1=1  i=l 

when  normalized  by  T  is  the  definition  of  the  TWAMP .   (TWAMP  as  defined  in 

Appendix  \    treats  time  as  discrete  in  units  of  months  and  therefore  must  make 

certain  adjustments  to  accomodate  installations  occuring  during  a  month.   The 
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differences  are  minor).   Thus,  Che  total  demand  over  the  program  time  base 
(0,T]  is  Poisson  distributed  with  parameter  aT*TWAMP.  Figure  4.1  illustrates 
this  for  the  case  in  which  9  total  installations  occur  over  the  interval 
(0,T]. 

If  all  of  the  installations  took  place  at  POC  (ti  =  t2  =  . . .  =  tjj  =  0) , 
the  process  describing  the  total  demand  is  a  homogeneous  Poisson  process  with 
rate  X    =  N  a.   However,  because  of  the  phased  installations,  the  total  demand 
process  has  a  mean  rate  which  is  a  nondecreasing  function  of  time;  i.e.,  the 

demand  process  is  a  non-homogeneous  Poisson  process.   The  mean  rate  of  demand 

t 
at  time  t  is  X(t)  ■  a  /   f(t)dt. 

0 
In   the   next    chapters  we   develop   alternative   models    for   determining    the 

amount    of    stock    to   be   provisioned   for   each   component    of    a   provisioning 

package.      In   all   of    those   models   we   assume   a   Poisson   distribution   for   the 

total   demand   for   each    item.      The   parameters    X±    referred    to   in   those   models 

will    be   assumed   to   be    the    "TWAMP'ed*'    values    discussed   above. 
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5.   UNITS  SHORT  MODELS  (SMA) 


INTRODUCTION 

Consider  a  weapon  system  composed  of  n   components.   Let   s.,s9,...,s 

be  the  number  of  units  of  components   1   through  n,   respectively,  to  be 

provisioned.   Let   c.   be  the  unit  cost  of  component   i   and  let   B  be  the 

total  budget  available  to  fund  spares  for  the  provisioning  package.   (B  will 

be  determined  by  the  COSDIF  process.)   As  discussed  earlier,  the  provisioning 

problem  is  interested  in  the  time  interval  (POC,  POC + TR+ PCLT) .   However, 

for  notational  convenience  we  will  refer  to  the  interval  as   (0,T.].   Let 

X.(T.)   be  a  random  variable  describing  the  total  number  of  demands  for  com- 
11 

ponent   i   in  the  interval   (0,T.);   p.(x)  =  Prob[X.(Tj  =  x] ;   P .  (x)  = 
Prob[X.(T.)  <_   x]  ;   and  H.(x)  =  Prob[X.(T.)  >  x]  .   As  indicated  in  the  previous 
chapter,  we  assume  that 

G.T.)Xexp(-A.T.) 

(    \  ii        ii 

p.(x)  =  j 

l  x! 

where   \.T.   is  the  total  expected  demand  over   (0,T.],   and   A.   is  estimated 

ii  ii 

from  the  technical  replacement  factor  (TRF)  and  the  installation  schedule. 


MODEL  I- -UN ITS  SHORT 

The  first  model  considers  the  objective  function  to  be  minimization  of  the 
expected  number  of  stockouts  over  the  interval   (0,Tj.   Thus,  we  want  to  find 
those  integer  values   s.,s?,...,s   which  minimize  the  total  expected  units 
short  subject  to  a  constraint  on  the  total  procurement  funds.   If   s.   units 
of  component   i   are  stocked  and  x.   units  are  demanded,  then  the  total 
number  of  units  short  for  conroonent   i   is 
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x  .  -  s  .     if   x  .  >  s  . 
11  11 


otherwise 


The  expected  number  of  units  short  for  component   i   is  therefore 


I  (x.  -  s  )  p  (x) , 

x  =s.+l    !ii 

1  1 


and  the  optimization  problem  is: 


n  °° 

(x 

l         ill 


Minimize   Z(S)  £  £  (x.-s.)    p.(x.),  (5.1) 


i=l   x.=s.+l 
i      l 


subject    to  2,      c.s.      £     B, 

i=l 


n 

c 
l    1 


and  s.  =  nonnegative  integer, 


This  constrained  nonlinear  separable  optimization  problem  can  be  solved 
using  the  generalized  Lagrange  multiplier  approach.   Let   S   be  the  n-vector 
(s.,s9,...,s  ),   and  let   L(S;9)   be  the  Lagrangian  function: 


n     °°  n 

L(S;8)   =  (x.  -  s.)p.(x.)  -  9(B  -     c.s.) 

•1         _i_1     !     Ill  .   ,    1  1 

1=1  X  =S .+1  1=1 

1    1 


Because  of  the  separability,  this  can  be  rewritten  as 


L(S;8)   =    I      L.(s.)  -  6B; 
i=l  x     1 


wnere 


L.(s.)   =  (x.  -  s.)p.(x.)  +  9c. s. 

11         L    ,  ,        1111      11 
x.=s .+1 
1   1 


If  the  demand  distributions  were  continuous,  the  optimal  solution  to  (5.1) 
would  be  given  by  solving  the  n+1  simultaneous  equations  for  s*,  i  =  1 
to  n,   and  0*  which  satisfy  [12]: 


3L. 

j^-     =      C  ;    i  =  1,2,  .  .  .  ,n, 
i 


^  =  0 

38 


Because  of  the  integer  nature  of  the  decision  variables,  the  above  procedure 
has  to  be  modified  slightly.   Let 


AL.(s.)   =  L.(s.)  -  L.(s.  -  1) 
li       li     li 


be  the  change  in  the  function  L.   resulting  from  an  increase  in  the  stockage 

level  for  component   i   from   s.-l   to   s..   The  optimal  value  of   s.   is 

11  l 

the  largest  value  of   s.   such  that   AL.(s)  <  0.   Simplifying,  we  determine 
that 


AL.(s.)   =   9c.  -  H.(s.) 

ii        ill 

Therefore,   s?   is  the  largest  nonnegative  integer  such  that 

H.(sp      :*c.     i  =  1,2, ...,n.  (5.2) 


The  above  simultaneous  equations  cannot  be  solved  in  closed  form  for  the 
decision  variables.   However,  the  solution  can  be  obtained  easily  using  the 
approach  by  Everett  (Reference  f  8  ]).   Select  a  trial  value  for  9   and  find 
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s*(6)   from  equation  (5.2)  for   i  =  l,2,...,n.   Then  evaluate 


n 

c 
i  i 


B(6)   =   I      c.s*(9) 
i=l 


If  B(9)  is  close  to  B  (within  prespecified  tolerance  limits)  the  process 
stops  and  s.  =  s . (9) .  If  B(9)  is  not  close  to  B,  a  new  value  of  9  is 
selected  and  the  process  repeats. 

Aid  in  the  selection  of  the  new  values  of  the  multiplier  is  available 
from  the  theory  of  the  generalized  Lagrange  multiplier  approach  presented  by 
Everett  [  8  ] .   The  multiplier   9   has  the  mathematical  interpretation  of  being 

9   "  "  SB  * 

Thus   9   can  be  thought  of  as  a  measure  of  the  decrease  in  the  total  units  short 
per  additional  budget  dollar.   Since  the  objective  function  is  a  monotone  non- 
increasing  function  of   B,   8   must  always  be  nonnegative.   If  the  amount  of 
a  current  solution  does  not  use  up  all  of   B   then  decrease   9   to  increase 
the  budget  expenditures;  if  it  exceeds  the  budget  then  increase   9   to  de- 
crease the  expenditures. 

The  theory  also  provides  guidance  at  each  iteration  as  to  the  quality  of 
the  incumbent  solution.   A  lower  bound  on  the  optimal  value  of  the  objective 
function  is  provided  (Daeschner  [5  ])  at  each  iteration  from  the  inequality 

Z(S*)   >   Z(S*(6.))  -  8.(B-B(8.)) 

J      J       J 

Comparison  of  the  present  value  of   Z;   namely,   Z(S*(9.)),   with  the  greatest 
lower  bound  yet  experienced  (for  cases  where  the  budget  constraint  was  satis- 
fied; i.e.,  the  amount  used  is   B(6.)    3)  will  provide  information  as  to 
the  maximum  benefit  to  be  gained  by  continued  searching.   For  example,  suppose 
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that  the  first  three  iterations  yielded  values  of  the  above  lower  bound  of 
(1003,  1155,  and  1167).   Also,  suppose  that  the  current  value  of  the  objec- 
tive function  is   Z(S*(9.))  =  1168.   We  know  that  the  optimal  solution  can- 

J 

not  be  smaller  than  1167,  so  there  is  little  to  be  gained  by  continued 
searching. 

The  iterative  process  described  above  will  determine  the  values  of  the 
decision  variables  which  will  minimize  the  total  expected  units  short  subject 
to  the  given  budget  constraint.   Exact  equality  of  budget  expenditures  may  be 
impossible  to  obtain  because  of  the  integer  nature  of  the  variables.   There- 
fore, the  result  determined  may  not  actually  be  optimal,  but  the  difference 
is  not  likely  to  be  significant.   It  is  important  to  appreciate  that,  for  a 
given  budget  expenditure  B(6),   the  solution  obtained  is  optimal. 

MODEL  II— MODIFIED  UNITS  SHORT 

The  model  of  the  previous  section  minimized  the  total  number  of  expected 
units  short  for  the  provisioned  system  subject  to  the  budget  constraint.   It  is 
a  trivial  matter  to  modify  that  model  to  allow  for  weighting  the  components 
by  essentialities.   Thus,  consider  the  objective  function 

n     °° 

Z0(S)   =  E.(x,  -  s.)  p.(x.) 

2        .u  u      ,   i   x    ill 

1=1  X .=s .+1 

1   1 

obtained  by  weighting  the  units  short  for  component   i   by  the  essentiality 
code   E.   where   E.  >  E.   if  component   i   is  judged  to  be  more  essential  than 
component   j.   The  minimization  of   Z_(S)   subject  to  the  budget  constraint 
gives  a  formula  for  finding  the  optimal  solution  which  is  similar  to  that 
found  in  the  previous  section.   The  optimal  stockage  level  for  item   i   is  the 
largest  nonnegative  integer   s.   such  that 


__ 


9c. 

H.(s.)   >   — - 
11       E. 


Similarly,  the  units  short  objective  function  can  be  modified  to  mini- 
mize the  total  expected  requisitions  short.   Let   g.   be  the  expected  quantity 
demanded  per  requisition  for  component   i,   and  consider  the  objective 
function 

n     a>    (x  .  -  s  .) 

Z3(S)   =    I  I  X    X   p.(x  )) 

i=l  x.-s.+l    Si 
1   l 

obtained  by  dividing  the  expected  units  short  by  the  expected  requisition 
size  to  yield  the  expected  requisitions  short.  For  this  problem,   s?"   is  the 
largest  nonnegative  integer  such  that 


H.(s.)   >   9c. g. 

i   i        i°i 


The  procedure  used  to  search  for  the  value  of   9   for  these  two  modifications 
is  the  same  as  that  discussed  for  Model  I. 

MODEL  III— SMA 

To  conclude  this  chapter  on  units-short  models  we  next  consider  an  objec- 
tive function  related  to  those  described  above,  but  which  is  stated  in  terms 
more  closely  related  to  presently  reported  measures  of  supply  effectiveness; 
namely,  supply  material  availability  (SMA). 

SMA  is  defined  as  the  fraction  of  all  requisitions  that  are  satisfied 
by  stock  on  hand  over  a  given  time  period.   It  would  seem  reasonable  to  at- 
tempt to  maximize  the  predicted  SMA  for  the  items  in  a  provisioning  package 
subject  to  a  constraint  on  the  total  procurement  investment. 

Let   D.   be  the  total  demand  for  item   i   during  the  provisioning 

interval   (0,Tj.   We  assume  that  requisitions  are  for  one  unit.   Let   F.   be 

i  1 


the  number  of  demands  satisfied,  and  let   U.   be  the  number  of  demands 

1 

unfilled.      Then,    clearly, 


D.      =      F.    +   U. 
1  11 


and  the  expected  number  of  demands  satisfied  (fills)  is 


E(F.)   =   E(D.)  -  E(U.)   =   a.T.  -  E(U.) 
l         l       l       ill 


Summing  over  the  n   items  in  the  provisioning  package,  we  get  the  total 
expected  number  of  "fills"  to  be 


n  n  n 

y  e(f.)   =    y  x.t.-  y  ecu.) 

. l-.  i        . l.  ii    .l,      i 

1=1  1=1  1=1 


n 

Division  by      X.T.   yields  the  ratio  of  expected  fills  to  expected  demands 

i=l 
( SMA) . 


n     <» 

V   C-/TT  >>  (x.  -  s.)p.(x.) 

E(F.)  L             L      .         l    l  ri   l 

>,     1  i=l  x.=s.+l 

SMA  =  — =1 — -1 (5.3) 

n  n 

)"   A.T.  y   a.T. 

.u ,      ii  >,   l  l 

i=l  i=l 


Observe  in  (5.3)  that  the  denominator  does  not  depend  on  the  decision  variables 
and  the  summation  term  in  the  numerator  in  the  total  expected  number  of  units 
short.  Thus,  the  allocation  of  a  given  budget  which  minimizes  expected  units 
short  will  be  the  same  as  the  allocation  which  maximizes  SMA.  However,  since 
the  value  of  SMA  is  likely  to  be  more  operationally  meaningful  than  the  total 
expected  units  short,  the  performance  prediction  should  probably  be  stated 
in  terms  of  SMA.   This  will  aid  in  the  communication  of  model  results. 
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Just  as  with  the  units  short  model,  the  SMA  model  can  be  modified  to 
include  essentiality  weights.   All  that  is  required  is  to  weight  the  units 
short  by  the  essentiality  values  and  to  modify  the  denominator  to  include 
also  the  essentiality  values.   The  essentiality-weighted  SMA  (ESMA)  would 
then  be: 


n  °° 

I      E             I           (x     -s    )p    (x) 
i=l      x   x.=s.+l        1111 
ESMA      =      1 - 


)  E  .  A  .  T  . 
.  **,  111 
i=l 


This  overall  value  can  also  be  written  in  terms  of  the  ESMA's  for  the  individual 
items.   Let  SMA.(s.)  be  the  expected  SMA  value  for  item  it  when  s.   units 
are  stocked.   Then  ESMA  can  be  written  as  follows: 

n 

y   E.A.T.  •  SMA.(s.) 
•=1   ill       11 

ESMA   =  — 

n 

J  E.A.T. 
.u,  ill 
i=l 


This  alternative  form  illustrates  that  the  overall  ESMA  value  is  a  weighted 
average  (weighted  by  essentialities  and  the  expected  number  of  demands,)  of 
the  individual  item  SMA's. 

SUMMARY 

In  this  chapter  we  have  considered  models  which  allocate  a  given  budget 
amount  so  as  to  optimize  system  performance  with  respect  to  units  short  or 
essentiality  weighted  units  short.   No  consideration  is  given  in  these  models 
to  the  length  of  time  that  must  pass  before  a  customer's  requisition  can  be 
satisfied  when  a  stockcut  occurs.   This  time  delay  is  uppermost  in  the  .ninds 


11 


of  the  customers  as  it  can  directly  affect  the  readiness  of  weapon  systems. 
The  models  developed  in  the  next  chapter  remedy  this  shortcoming  by  explicitly 
considering  time-weighed  units  short. 
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6.   TIME-WEIGHTED  UNITS  SHORT  MODELS 

We  develop  two  models  in  this  chapter  which  consider  not  only  the 
shortages  which  accumulate  over   (0,T],   but  also  the  length  of  time  that 
each  shortage  exists.   Classical  inventory  cost  models  recognize  the  effect 
that  stockout  delays  have  on  system  performance  by  imposing  a  stiff  penalty 
cost  for  stockouts  which  grows  as  the  stockout  time  increases.   The  models 
that  we  develop  here  do  not  assess  penalty  costs  since  they  are  virtually 
impossible  to  quantify.   Instead,  the  objective  function  is  the  direct  mini- 
mization of  expected  time-weighted  units  short  (TWUS) .   Such  an  objective 
function  is  appropriate  when  the  readiness  of  the  operational  units  suffers 
as  long  as  replacement  parts  are  unavailable.   In  a  readiness  measure  such 
as  operational  availability,  the  time  until  a  failed  component  is  restored 
to  serviceable  condition  is  a  key  ingredient.   When  spares  are  available,  that 
time  is  short.   But  when  spares  are  not  available,  that  time  will  be  very 
long. 

These  models  assume  that  if  there  is  a  spare  item  in  the  system  then 
it  is  "immediately  available"  and  that  if  there  is  no  spare  the  demand  be- 
comes a  backorder  which  is  filled  by  the  next  order  which  arrives  from  the 
manufacturer . 

MODEL  IV— TWUS 

Let  us  restrict  attention  initially  to  a  single  component.   Let  X(T) 
be  the  random  number  of  demands  for  the  component  in   (0,Tj .   Let  us  assume 
that   X(T)  =  m,   and  let 

0  <  T..  <  T_  <  .  .  .  <  T   <T 

12  m  — 

be  the  times  at  which  the  ra  demands  occur.   If   s   units  of  cue  icem  are 


stocked,  the  total  time-weighted  units  short,  given  X(T)  =  m  and 

T, ,T„,...,T  ,   can  be  described  by 
1  2      m 


m-s 
[   (T-T  ,  .). 

3-1     S+J 


Figure  6.1  provides  a  graphic  illustration  of  this  derivation.   The  backorders 
are  those  demands  giving  negative  net  inventory.   The  length  that  each  is 
outstanding  is  the  time  from  when  its  demand  occurs  until  the  next  order 
arrives  at  time  T. 


NET 
INVENTORY 


m-s 


Figure  6.1:   NET  INVENTORY  OVER  TIME 


The  conditional  expected  value  of  the  time-weighted  units  short  is  therefore: 


if  m  <  s 


E[TWUS(s)  JX(T)  =m] 


(6.1) 


m-s 


(T-E[T    X(T)  =m])   if   m  >  s 
j=l       S+J 


2S 


Let  us  now  determine  E[T,  |X(T)  = m] .   Assuming  a  homogeneous  Poisson  process, 
the  probability  that  the  kth  demand  will  occur  between  t   and   t+dt   can 
be  determined  as  follows: 


f_  (tlm)dt   =  Prob(t  <  T,  <  t+dt|x(T)  =  m) 
T,  K 


k 


P(k-1  events  in  (0,t],  1  event  in  (t,t+dt],  m-k  events  in  (t,T]) 

P(X(T)  =  m) 

(At)k"1exp(-At)n    ,(A(T-t))m"kexp(-A(T-t)) 
(k-1)!      (A  dt)        (m-k)! 

(AT)mexp(-AT) 


m!        /t.k-1  ,T-t,m-k  dt 


(k-1) !  (m-k) !  VT'     v  T  '     T  ' 

Now,   E[T,  |X(T)  =  mj  =  jl    t  f_  (t|m)dt 
k  '  (J     T, 


J0   C  (k-1)!  (m-k) !  V     "   T    }  T 


The  change  of  variable,   u  =  t/T,   simplifies  the  above  expression  to 


_  (1  m!         k-1  ,,    m-K. 

T  lA  u  -r. — 7— t— 7 r^-r  u     (1-u)     du 

;  0    (k-1)  !  (m-tc)  ! 


The  integral  is  recognized  to  be  the  expected  value  of  a  beta  distributed 
random  variable  with  parameters  a  =  k  and  3  =  m-k+1.   Consequently,  the 
conditional  expectation  that  we  seek  is 


T  (-?*■> 


'a+6      m+1 


Returning  co  expression  (6.1),  we  have  determined  the  conditional  expectation 


m-s 

T  _  T  I 

^m+1 
E[TWUS  X(T)  =m]   = 


I 


I      T-T(gf)   if   m  >  s 


if  m  <  s 


Expanding  the  summation  gives  (for  the  case  m  >  s) 

T(m-s) (m-s+1) 
2(m+l) 

Finally,  on  summing  over  all  possible  values  of  m,   the  total  expected  time- 
weighted  units  short  for  the  component  is  found  to  be 


E[TWUS(s)]   =    )   E[TWUS(s)  |N(T)  =  m]  Prob[N(T)=m] 
m=0 


co  in 

T  (m-s)(m+l-s)  (AT)  exp(-XT) 

'■  ,  2     nri-1  m! 

m=s+l 


|{(AT)H(s)  -  2sH(s+l)  +  S(S^1)  H(s+2)} 


|{H(s+l)[AT-2s  +  s(^1)]  +  p(s)(AT-s)} 


Now  that  we  have  derived  an  expression  for  the  expected  time-weighted 
units  short  as  a  function  of  the  stockage  level   s,   let  us  write  the  optimi- 
zation problem  that  we  want  to  solve.   Let   V.(s.)   be  the  expected  time- 
weighted  units  short  for   item   i  when  s.   units  are  provisioned.   We 
want  to  minimize  the  total  essentiality  weighted,  time-weighted  units  short 
subject  to  a  budget  constraint.   The  Lagrangian  statement  of  the  problem  is: 


n  n 

Minimize   L(S;6)   =   T   E . V . (s . )  -  G [3  -   /   c.s.] 

i=l  L   X  X  i=l   X  ± 
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Let   L  (s  ;9)  =  E.V. (s.)  +  3c. s..   As  before,  the  optimal  solution 
11       111      11  '      r 

s^   is  found  by  determining  the  largest   s.   such  that   AL.(s.;9)  <  0  wher< 

I      c  s.  =  B.   After  simplification  AL.(s.;6)   can  be 
i=l  x     x 


6   is  such  that 
written  as: 


a  . 

AL.(s.;9)   =  Ei[(l-P(si))(Ti  -  ~)    +  T.p(s.)]  -  8c. 


Thus,      s.      is   the    Largest      s_.      such   that 


s  .  0c  . 

(l-P.(s.))(T,    -  t^)  +  T.p.(s.)    <  -rA 

ll           l         A.  ill             E. 

l  l 


We  solve  for  the  optimal  values  using  exactly  the  same  search  (on  3) 
procedure  as  discussed  in  the  previous  chapter. 

MODEL  V— MSRT 

One  of  the  most  widely  reported  measures  of  supply  performance  is  the 
mean  supply  response  time  (MSRT) .   This  measure  inherits  its  popularity  at 
least  partly  through  the  role  it  plays  in  the  expression  commonly  used  to 
define  availability  (see  Chapter  7).   In  addition,  MSRT  is  important  by  it- 
self as  an  indicator  of  the  success  of  the  supply  system  in  meeting  response 
time  goals.   As  a  measure,  it  considers  both  the  likelihood  of  satisfying 
demands  from  stock,  on  hand  and  the  lengtn  of  the  delay  in  satisfying  demands 
wnen  the  system  runs  out  of  stock. 

Let  us  determine  the  allocation  of  a  fixed  budget  which  will  minimize 

the  MSRT  averaged  over  all  of  the  items  in  a  provisioning  package.   Let   t 

be  the  time  at  which  a  failure  occurs  in   (G,T.j,   and  let   Y.(t)   be  the 

total  number  of  failures  which  have  occurred  from  FOC  up  to,  but  not  ^ncluaing, 

t.   Let   R.(t)   be  the  supply  response  time  for  the  aemand  occurring  at  t±rr,e 
i 

t.   We  assume  that   R-(t)   is  a  constant,   k.   (possib.lv  zero)  if   Y.(t)  <  s.- 

L  1  11 
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If  the  on-hand  stock  is  not  positive  when  the  demand  occurs  at  time   t,   then 

satisfaction  of  the  demand  will  be  delayed  until   k.   units  of  time  bevond 

1 

receipt  of  the  first  replenishment  order  which  we  assume  will  occur  at  time 
T.   Thus, 


R±(t)   = 


k.         if  Y. (t)  <  s. 
1  11 


k.  +  (T.-t)   if  Y  (t)  >  s. 

li         i    —  l 


The  conditional  mean  supply  response  given  the  failure  occurs  at  time   t  is 


E[R.(t)]   =  k.P[Y. (t)  <  s. J  +  (k.  +  (T-t))P[Y.(t)  >  s. 
l  li      l      li       l—i 


=  k.  +  (T-t)P[Y.  (t)  >  s.] 
l     -;      1-1 


The  unconditional  mean  supply  response  time  is  found  by  integrating  the  above 
conditional  expectation  over   (0,T]   with  respect  to  the  density  of  the  time 
t.   But,  it  is  well  known  (see  Ross  [21])  that  the  time  of  failure  of  an 
arbitrary  event  generated  by  a  Poisson  process  is  uniformly  distributed  on 
the  interval   (0,T].   Thus,  the  density  of   t   is 

ft(T)   -      - 

'  0   otherwise 

and  the  mean  supply  response  time  is 

T. 

MSRT.(s.)   =   f   (k.  +  (T-t.)H.(s.))  -  dt 
11      '  i     1111    T. 

0  i 

T. 
i 


=   k.  +  \  j      (T-t.)H.(s.)dt 
l    T.  ;  .  x     i   i   i 


=   k 


i  0 

T.  s.  (s.+l)s. 

. +-k.(s.)  -  —  H.(s.+1)  +  — ^ ^  H.(s.+2) 

i  2  i  i  A.ii  <->  -      i  x 


A.  T. 
l  i 
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Recalling  the  expression  for  the  expected  time-weighted  units  short  we 
observe  that 

V  (s  ) 

MSRT.(s.)   =  k.  + 


i   i       i     A.T. 
1  1 

where  V  (s  )   is  the  expected  time-weighted  units  short  when  s    units  of 
11  l 

item   i   are  stocked.   Finally,  the  aggregate  mean  supply  response  time  for 
all  items  in  the  provisioning  package  is  obtained  by  weighting  the  individual 
means  by  the  expected  cumulative  demand:  * 

n 

y  (a.t;msrt.(s.) 


MSRT   = 


n 

I      A.T  . 
i=l  X    X 


n  n 

I      k.A,T.+     V.(s.) 
.  ,   ill   ,t1   l   l 
1=1 i=l 

n 

I      A.T. 
i=l   X  X 


(6.2) 


Since  the  denominator  and  the  first  summation  in  the  numerator  do  not 

depend  on  the  decision  variables,   sn   through   s  ,   it  is  clear  that  the 

1       c  n 

stockage  levels  which  minimize  the  aggregate  MSRT  are  the  same  as  those  that 

minimize  the  total  expected  time-weighted  units  short.   As  with  the  TWUS  model, 

the  MSRT  model  can  be  modified  easily  to  incorporate  essentiality  weights. 

All  that  is  needed  is  to  replace  the  weights   A.T.  in  (6.2)  by   E.A.T.  . 

r  °  li  'ill 

Even  though  the  TWUS  model  and  the  MSRT  model  will  yield  the  same  solu- 
tions, MSRT  appears  to  be  more  operationally  meaningful  than  is  time-weighted 
units  short.   The  next  chapter  develops  a  model  which  incorporates  MSRT  to 
yield  another  measure  of  performance  of  interest  to  the  operational  community; 
namely,  operational  availability. 
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7.   AVAILABILITY  MODEL 

The  final  model  which  we  develop  as  a  candidate  for  provisioning  is  an 
availability  model  which  incorporates  as  a  key  ingredient  the  MSRT  discussed 
in  the  previous  section.   We  believe  that  the  Navy  Supply  System  would  be 
better  served  by  using  one  of  the  previous  models,  but  we  include  the  availa- 
bility model  here  to  satisfy  the  present  clamor  for  models  which  "spare  to 
availability". 

The  commonly  used  definition  of  the  availability  for  a  single  component 
is 

A.   =   MTBF./(MTBF.  +MTTR.  +MSRT.) 
l         1111 

where  MTBF .   is  the  mean  time  between  failures,   MTTR.   is  the  mean  time  re- 
l  i 

quired  to  repair/replace  a  failed  unit  when  a  spare  is  available,  and  MSRT. 
is  the  mean  supply  response  time  {see   References  [  3  ],  [14],  [15],  [23]). 
Although  this  expression  is  actually  a  limiting  result  which  is  valid  only 
under  certain  conditions  ([13],  [21]),  it  is  useful  for  providing  insights 
into  the  factors  which  affect  availability  and  of  the  sensitivity  of  availa- 
bility to  those  factors.   For  purposes  of  provisioning,  we  write  the  mean 
supply  response  time  as  a  function  of  the  total  number  of  units  provisioned, 
MSRT.(s.),   as  was  done  in  the  previous  section.   The  availability  expression 
then  becomes 


A.(s.)      =      MTBF./ (MTBF.    +MTTR.    +MSRT.(s.))  (7.1) 

li  li  l  li 


Clearly,  a  primary  objective  of  the  Supply  System  should  be  to  allocate 
resources  so  that  weapon  system  availability  is  maximized.   The  key  phrase 
^n  the  previous  sentence  is  "weapon  system."  This  is  where  one  of  the  major 
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difficulties  arises  in  trying  to  develop  availability  models   The  availability 
of  a  weapon  system  is  a  very  complicated  function  which  depends  on  1)  the 
availabilities  of  its  components,  2)  the  configuration  of  the  weapon  system 
(reliability  block  diagram),  and  3)  the  way  that  the  system  is  deployed 
operationally.   Without  making  some  simplifying  assumptions,  it  would  probably 
not  be  possible  to  develop  any  useful  general  purpose  algorithms  for  maxi- 
mizing weapon  system  availability. 

Most  existing  availability  models  (see  References  [  3  ]  and  [23] ,  for 
example)  make  the  following  assumptions: 

1)  Component  failures  are  generated  by  a  Poisson  process  (lifetimes 

are  exponential).   This  means  that  MTBF .  =  1/A . . 

i      l 

2)  Component  availabilities  are  given  by  expression  (7.1). 

3)  The  weapon  system  is  configured  with  the  components  all  in  series. 

4)  The  failure  of  one  component  is  independent  of  the  failures  of 
the  others. 

These  assumptions  lead  to  Che  formula  below  for  the  availability  of  the 
weapon  system. 


A  (s   s0,...,s  )   =    H  A.(s.)  (7.2) 

o   1   2      n       .,ii 

1=1 


The  accuracy  of  the  above  expression  depends  on  the  extent  to  which  the 
assumptions  are  valid.   We  doubt  that  they  are  all  satisfied  for  any  real 
weapon  system,  and  we  doubt  that  the  resulting  estimated  weapon  system  availa- 
bility would  accurately  reflect  the  actual  observed  system  availability.   As 
a  consequence,  we  will  refer  to  this  model  as  the  "pseudo- availability  model." 
Nevertheless,  as  a  model  of  system  availability,  (7.2)  may  be  useful  for 
the  purpose  of  deciding  the  range  and  depth  of  items  to  be  provisioned. 
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Consider  the  problem: 


n 

maximize     A        =        II     MTBF./(MTBF.    +  NTTTR.    +MSRT.(s.))  (7.3) 

o      .  ,      i       i        1        11 
sr...,sn        1=1 


n 

s.t.      J   c.s.   <   B  s.  =  0,1,2,  .  .  . 

.**,   l  l  —  l 

1=1 

Taking  the  logarithm  of  both  sides  of  the  objective  function,  we  get: 


n 

inA(8,,s, )   =      £n(MTBF./(MT3F.  +  MTTR.  +  MSRT . (s . ) ) ) 

1     £.  n      .Ln        l      i       i       11 

i=l 


Let      L(S.;9)    =   £n(MTBF  ./MTBF.  +ffTTR.  +MSRT  .  (s  . ) )    -    9c.  s..      Solving   problem 
l  11  l  11  11  &r 

(7.3)    is   eauivalent    to   solving   the  Lagrange   problem: 


n 

max   L(s,,...,s    ;0)      =        V     L.(s.;9)   +  OB. 
1  n'  .*•_      li 

i=l 


The  optimal  solutions  are  given  by  determining  the  largest  values  of   s. 
such  that 


AL.(s.;9)   =  L.(s.;9)  -  L.(s.-1;9)   >   0. 
11         11       11 

Simplification   shows    that      s?      is    the   largest    non-negative   integer   such   that 

J.    +  MSRT.(s.-l) 

J1  +  MSRT\(s  J  exp(9*c.), 

l  11 

where     J.    =  MTBF.+MTTR.      and      9*      is    that    nonnegative   Lagrange  multiplier 
ill  ° 

n 

value  such  that    )   c.s?  =  B. 

1-1   X  L 

Unlike  the  previous  models,  it  would  not  make  sense  to  weight  the  indi- 
vidual components  by  their  essentialities  in  the  availability  model.   In 
fact,  the.  series-system  assumption  is  tantamount  to  an  assumption  of  equal 
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essentiality  of  components;  i.e.,  the  failure  of  any  component  will  render 
the  system  unavailable.   True  availability  models  which  explicitly  consider 
system  configuration  automatically  handle  essentiality  by  consideration  of 
the  impact  of  each  component  failure  on  system  performance  (see  Jee  [13]). 
No  external  weighting  or  assignment  of  essentiality  cedes  would  be  necessary. 
Such  "true"  availability  models  are,  however,  very  difficult  to  develop  and 
implement . 

Because  of  Che  assumptions  made  in  developing  the  pseudo-availability 
model  in  this  section,  it  is  unlikely  that  the  system  availability  predic- 
tions provided  by  the  model  will  be  correct.   Indeed,  for  even  relatively 
small  systems,  the  predicted  availability  will  necessarily  be  very  small. 

For  example,  a  system  composed  of  25  components,  each  with  an  availability 

25 

of  0.95,  will  have  a  predicted  availability  of  only  (.95)"   =  0.2774. 

Nevertheless,  the  budget  allocation  decisions  made  by  the  model  may  be  very 
good. 

To  conclude  this  discussion  about  the  pseudo-availability  model  we  would 
like  to  comment  on  a  claim  frequently  made  in  the  availability  literature 
(see,  for  example,  [  3]  and  [23]).   Some  of  the  availability  analysts  argue 
that  when  sparing  to  availability,  the  only  term  in  the  availability  expres- 
sion which  involves  the  number  of  spares  is  MSRT.   Therefore,  they  claim, 
maximization  of  system  availability  with  respect  to  the  number  of  spares 
for  the  components  is  equivalent  to  minimization  of  MSRT.   A  simple  counter- 
example given  in  Appendix  C  shows  that  the  claim  is  not  true!   However,  for 
realistic  choices  of  the  various  parameters,  the  allocacions  determined 

from  the  two  models  CMSRT  and  A  )  may  be  very  much  alike.   The  numerical 

o 

evaluations  which  follow  in  the  second  reoort  examine  this  issue. 
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8.   COMMENTS  ON  THE  MODELS 

INTRODUCTION 

The  previous  chapters  have  presented  three  basically  different  provision- 
ing models: 

(1)  units  short  (SMA) , 

(2)  time-weighted  units  short  (MSRT) ,  and 

(3)  pseudo-availability. 

In  this  chapter  we  provide  brief  comments  about  the  models  as  related 
to  three  areas: 

(1)  incorporation  of  a  random  protection  interval, 

(2)  theoretical  properties  of  the  models,  and 

(3)  solution  of  a  related  problem  for  use  as  a  budgeting  tool. 

RANDOM  PROTECTION  INTERVAL 

The  previous  models  all  assume  that  the  length  of  the  protection  inter- 
val (POC  to  POC+T]  is  known.   In  reality,   T   is  a  random  variable  since  both 
the  procurement  leadtime  and  the  time  from  POC  to  the  first  replenishment  buy 
are  random  variables.   If  one  had  information  about  the  distribution  of  the 
random  variable  T,   the  previous  models  could  be  extended  (at  the  cost  of 
additional  mathematical  complexity)  to  accommodate  the  randomness  of  T. 

Let   gT(-)   be  the  probability  density  function  for  T.   What  is  required 
in  the  previous  models  to  incorporate  the  randomness  of  T   is  to  replace  T 
everywhere  by  a  dummy  variable,  say   u,   and  to  integrate,  weighting  by  the 
density   g^Cu),   over  all  possible  values  of   u.   One  common  assumption  for 
the  form  cf   grp(u)   is  the  gamma  density.   The  gamma  probability  family  has 
the  proper  range  of  values   (0   to  °°) ,   and  the  family  is  extremely  rich 
in  the  shapes  that  it  can  assume.   Many  theoretical  and  empirical  studies 
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have  shown  that  the  gamma  distribution  frequently  approximates  well  the 
empirical  distributions  of  leadtime.   Furthermore,  when  the  gamma  distribu- 
tion for  leadtime  is  combined  with  the  Poisson  distribution  for  demand  in  an 
interval,  the  resulting  distribution  of  leadtime  demand  is  negative  binomial 
(see  Chapter  3  of  [12]).   This  theoretical  fact  accounts  for  the  selection 
of  the  negative  binomial  as  one  of  three  distributions  used  by  the  Navy  ICP's 
in  some  of  its  wholesale  provisioning  and  replenishment  models. 

For  purposes  of  implementation  in  the  provisioning  models  it  is  unlikely 
that  any  information  will  be  available  to  adequately  estimate  the  distribution 
of  T.   Also,  since  inclusion  of  a  probability  distribution  for  T   increases 
the  complexity  of  the  models,  we  suggest  that   T   continue  to  be  treated  as 
a  constant.   Of  course,  every  effort  should  be  made  to  estimate  T   as  well 
as  possible. 

One  of  the  elements  of   T   is  the  length  of  time  from  POC  to  the  time 
that  the  first  replenishment  buy  is  made.   In  addition  to  the  elements  already 
included  in  the  models,  this  time  period  depends  on  an  external  element  not 
yet  considered — the  replenishment  reorder  level.   Given  the  size  of  the  first 
provisioning  buy,  say   s,   and  the  size  of  the  first  replenishment  reorder 
level   r,   it  is  easy  to  determine  that  the  time  from  POC  to  the  first  replenish- 
ment buy  is  Erlang  distributed  (assuming  a  homogeneous  Poisson  process  with 
rate   A   generates  demands)  with  parameters   n  =  s-r   and   \.   This  is  because 
the  delay  is  the  amount  of  time  required  to  experience   s-r   demands  with 
exponentially  distributed  times  between  demands.   The  problem  is  chat  in  order 
to  integrate  the  provisioning  and  replenishment  models  it  would  be  necessary 
to  include  the  replenishment  reorder  level   r   as  yet  another  variable  in  the 
provisioning  model  so  that  the  size  of  the  provisioning  buy  depends  on   r. 
Furthermore,  the  value  of   r   in  such  an  incegrated  model  would  logically 
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depend  on  the  provisioning  stockage  level  so  that  solution  of  the  model  would 
require  an  iterative  technique  greatly  complicating  the  model  and  increasing 
the  computational  burden.   Thus,  we  suggest,  that  the  first  replenishment 
buy  not  be  triggered  by  the  inventory  position  dropping  below  a  reorder 
point,  but  be  made  automatically  at  a  fixed  time  after  POC,  say  POC  +  At. 
If  the  estimates  of  demand,  installation  schedules,  and  procurement  leadtimes 
are  accurate,  it  is  reasonable  that  the  first  replenishment  buy  be  made  at  POC 
(i.e.,   at  =  0) .   Delaying  the  replenishment  buy  to  some  later  date  to  allow 
the  initial  TRF  estimates  to  be  updated  by  actual  demand  data  only  lengthens 
the  provisioning  protection  interval  and  increases  the  size  of  the  required 
provisioning  buy. 

THEORETICAL  PROPERTIES 

The  models  developed  in  this  report  can  be  expected  to  lead  to  different 
solutions  to  the  optimal  budget  allocation  problem.   In  contrast  to  the  Varia- 
ble Threshold  Model,  each  of  the  models  also  provides  a  prediction  of  the 
expected  performance  to  be  achieved  from  a  given  solution.   Furthermore,  some 
theoretical  properties  of  the  models  can  be  examined  to  determine  which  model 
best  reflects  our  intuitive  notion  about  how  the  system  is  degraded  by  the 
presence  of  shortages,  about  the  effect  of  uncertainty  in  our  estimates  of 
the  demand  rates,  and  about  the  propensity  to  stock  items  as  a  function  of 
their  demand  rates,  unit  costs,  essentialities,  and  leadtimes. 

All  of  the  models  are  more  likely  to  stock  the  low  cost,  high  demand-rate 
items,  all  other  things  the  same.   Weighting  by  essentialities  is  the  only 
way  offered  here  to  oppose  that  tendency.   Also,  the  models  are  all  piece- 
parts  oriented.   No  explicit  consideration  is  given  in  any  of  the  models  for 
"weapon  system  effectiveness".   However,  essentiality  weights  can  be  used  to 
reflect  the  importance  of  the  individual  items  with  respect  to  weapon  system 
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effectiveness.   Proper  assignment  of  essentiality  weights  should  reflect 
what  happens  to  the  readiness  of  a  weapon  system  if  a  shortage  occurs  in  an 
individual  item. 

A  COMPANION  PROBLEM  AS  A  BUDGETING  TOOL 

In  the  models  which  we  have  developed,  we  have  assumed  that  a  budget   B 
was  given  at  the  outset.   DODI4140.42  requires  that  the  provisioning  budget 
be  determined  by  using  the  COSDIF  model.   As  discussed  in  Appendix  B,  the 
COSDIF  model  first  determines  if  an  item  is  to  be  stocked  by  comparing  the 
expected  cost  of  stocking  the  item  with  the  expected  cost  of  not  stocking  the 
item.   If  the  former  cost  is  smaller,  the  item  satisfies  the  range  rule. 
The  depth  of  a  selected  item  is  the  expected  demand  during  a  procurement 
leadtime  plus  one  quarter.   The  total  "priced-out"  value  for  all  items  satis- 
fying the  range  rule  is  the  provisioning  budget  (excluding  the  budget  allowed 
for  insurance  and  NSO  items) . 

This  procedure  for  determining  the  budget  for  a  provisioning  budget  is 
cost-based  requiring  cost  estimates  that  may  not  reflect  accurately  real 
costs.   Furthermore,  the  COSDIF  procedure  does  not  provide  any  estimate  of 
the  performance  that  might  be  expected  from  the  resulting  budget.   The. 
alternative  models  presented  in  this  report  can  each  be  rewritten  in  a 
mathematically-related  form  which  could  be  solved  to  provide  a  tool  useful 
for  determining  a  provisioning  budget. 

We  saw  that  each  of  the  models  described  earlier  could  be  written  in 
the  primal  form: 

n 

min(max)    \      f  .  ( s  . ) 

1   i. 


i=l 


lbiec:  to      >   c .  s  .   <   B 

,  -  ,   i  i  — 


(8.1) 


The  companion  problem  can  be  defined  as  follows: 


minimize    )  c.s. 
11 


subject  to       f.(s.)  <   (_>)   G' 
i-1  ±      X 


where  G*   is  a  specified  performance  goal.   (The  direction  of  the  inequality 
depends  on  whether  (8.1)  is  a  minimization  or  a  maximization  problem).   This 
companion  problem  can  either  be  solved  directly  using  a  generalized  Lagrange 
approach  like  that  used  to  solve  the  primal  problems,  or  the  primal  problem 
can  be  solved  for  a  range  of  budget  values  and  the  smallest  budget  satisfying 
the  performance  goal  be  selected. 

Solution  of  this  companion  problem  will  provide  planners  with  a  budgeting 
tool  through  which  they  can  justify  a  budget  on  the  basis  of  the  performance 
that  can  be  attained.   Thus  it  provides  a  mechanism  for  relating  resources 
to  readiness. 


9 .   SUMMARY 

We  have  developed  alternative  wholesale  provisioning  models  which  are 
intended  for  use  by  both  SPCC  and  ASO.   The  models  determine  how  to  allocate 
optimally  a  given  provisioning  budget  with  respect  to  various  measures  of 
performance.   The  measures  of  effectiveness  can  be  summarized  into  three  types: 

1)  Essentiality-weighted  units/requisitions  short  or  SMA, 

2)  Essentiality-weighted  time-weighted  units  short  or  MSRT,  and 

3)  Availability. 

In  addition  to  providing  the  optimal  budget  allocations,  each  model  provides 
a  prediction  of  the  performance  to  be  expected  during  the  provisioning  period. 

The  models  all  assume  that  a  budget,  as  determined  by  the  COSDIF  process, 
is  given.   However,  each  model  can  be  converted  easily  into  a  related  optimiza- 
tion problem  which  determines  the  minimum  budget  required  to  achieve  a  speci- 
fied performance  goal.   We  feel  that  it  makes  more  sense  to  determine  a  budget 
in  this  manner,  than  by  using  the  COSDIF  approach  which  is  heavily  cost 
oriented  and  which  provides  no  estimates  of  the  performance  to  be  expected 
with  a  given  budget. 

The  alternative  models  solve  simultaneously  the  range  and  the  depth 
problems.   Thus,  no  additional  model  is  required  for  the  budget  allocation 
process . 

In  order  to  evaluate  the  proposed  alternative  models,  the  modeis  should 
be  exercised  with  sample  data,  and  the  results  examined  carefully  to  observe 
how  each  model  allocates  the  provisioning  budget.   The  results  of  each  model 
should  be  tested  with  respect  to  several  measures  of  performance  (units  short, 
SMA,  time-weighted  units  short,  MSRT,  availability,  etc.;.   One  should  then 
select  that  model  which  provides  results  that  best  agree  with  a  decision 
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maker's  perception  of  how  the  system  should  perform.   In  exercising  the 
alternative  models  the  present  models  should  also  be  tested  to  serve  as  a 
baseline  for  comparison.   This  is  the  type  of  model  evaluation  that  will  be 
discussed  in  the  second  report. 

However,  even  before  a  real  data  performance  evaluation  is  performed, 
the  alternative  models  can  be  evaluated  theoretically  or  intuitively  to  deter- 
mine if  they  attempt  to  achieve  the  goals  of  the  supply  system.   We  have  com- 
mented on  various  properties  of  each  model  in  the  sections  developing  the 
models  and  in  the  previous  chapter.   Because  of  those  theoretical  properties 
of  the  models,  and  because  of  the  role  that  time-weighted  shortages  (MSRT) 
plays  in  the  expression  for  availability,  we  believe  that  the  MSRT  model  makes 
the  most  sense  for  provisioning.   It  captures  the  essence  of  the  objectives 
of  the  supply  system  while  avoiding  some  of  the  problems  and  data  requirements 
that  would  be  associated  with  the  pseudo-availability  model. 
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APPENDIX  A:   INITIAL  DEMAND  FORECASTING 

THE  PRESCRIBED  PROCEDURE 

The  initial  demand  prediction  process  has  been  delineated  in  Reference 
[7]  and  begins  with  a  forecast  of  "program  data"  from  POC  out  over  the  two- 
year  demand  development  period.   This  data  is  a  month-by-month  schedule  of 
the  anticipated  end  item  installations  over  that  time  period.   The  time- 
weighted  average  months  program  (TWAMP)  is  computed  next  for  the  appropriate 
program  time  base  (PTB) .   Reference  [7]  allows  the  PTB  to  have  a  maximum 
length  of  one  year.   Shorter  PTBs  are  3  months  and  6  months  and  are  appro- 
priate for  computing  needs  for  expensive  spares.   However,  SPCC  and  ASO  use 
only  a  12  months  PTB  regardless  of  an  item's  value  because  of  the  workload 
created  by  having  shorter  PTBs  (a  review  is  needed  at  the  end  of  each  PTB 
interval  for  the  shorter  ones) . 

The  concept  behind  the  TWAMP  is  quite  simple.   Consider  Figure  A.l;  the 
TWAMP  is  computed  by  determining  the  area  under  the  curve  of  the  installation 
schedule  and  chen  dividing  it  by  the  length  of  the  PTB.   This  results  in  an 
average  number  of  end  items  to  be  supported  over  the  PTB.   This  average  is 
denoted  as  the  initial  TWAMP  or  TWAMP  by  Reference  [19]. 


Fi"ure  A.I:   An  Installation  Schedule 
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The  following  formulas  are  provided  by  Reference  [7]  for  computing  the 
initial  TWAMP .   They  obtain  the  area  by  summing  the  vertical  slices  of 
the  area  spanning  each  month.   Any  installations  occurring  in  a  specific 
month  are  assumed  to  begin  operation  in  the  middle  of  the  month. 


D 
m 


11 

-r-  for   m  =  1 


m-1      I 

I      I  +  -f        for  m  >  2 
k=l   k 


where 


I   =   the  number  of  units  of  the  end  item  to  be  installed 
in  month  k. 


PTB 

p 

I 

J 

m=l 

m 

TWAMPj       pTB   . 

The  initial  annual  demand  rate  for  a  given  spare  part  can  be  computed 
by  using  the  following  formula  [19]: 


D,   =   TWAMP _  *  N  *  3RF 


where 


N   =   the  number  of  units  of  a  given  replaceable  part  in 
the  end  item. 

5RF  =   tne  "best  replacement  factor"  which  is,  in  fact,  the 
estimated  failure  rate  of  one  unit  over  a  year. 


The  BRF  is  initially  based  on  a  technical  replacement  factor  (TRF)  which 
is  the  contractor's  estimate  of  the  attrition  rate. 

The  COSDIF  formula  requires  an  initial  estimate  of  the  steady  state 
annual  demand  rate  D,^.   This  is  computed  by  first  determining  the  number 
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of  end  items  to  be  installed  over  the  first  two  years  (its  value  is  denoted 
as  TWAMP    in  Reference  [19]).   Then 

Dss   =   TWAMP   *N*BRF. 

If  a  spare  part  is  repairable  then  the  formula  considers  only  attrition 
demand.   In  that  case, 


Dgs   =   TWAMPSS*N*BRF[1  -  CRR*RSR]  (A.l) 

where: 

CRR  =   the  carcass  return  rate,  and 
RSR  =   the  repair  survival  rate. 

The  initial  estimate  of  demand  for  repairables  will  be  introduced 
later. 

The  development  of  the  provisioning  budget  constraint  of  Reference  [7] 
requires  that  the  initial  demand  during  the  procurement  leadtime  (PCLT) 
plus  one  quarter  be  computed.   This  is  then  the  depth  value  to  be  used  to 
develop  the  budget  constraint.   Reference  [19]  also  needs  the  demand  during 
PCLT  for  the  Variable  Threshold  model.   For  consumable  items  the  formulas 
are,  respectively, 


D       =   D  M*; 

PCLT      A   4  ' 


=    ,(PCLT+1) 

PCLT+1      A     4 


For  a  repairable  item  the  computation  used  in  Reference  [19]  includes  the 
repair  turnaround  time  (TAT).   The  formulas  are  based  on  the  following  logic. 
Without  any  carcass  returns  the  attrition  demand  is  the  same  as  Equation 


k  : 


(A.l)  with  TWAMP  „<,   replaced  by  TWAMP  .   If  the  repair  process  is  "up  and 
running"  at  POC  minus  TAT,   then  the  first  repaired  item  will  be  available 
at   POC  and  will  provide  a  replacement  for  some  unit  of  the  item  when  it 
fails.   The  net  attrition  loss  is  then 

TWAMP  *N*BRF*[1  -  CRR*RSR] 

If,  as  is  expected,  the  repair  process  can't  begin  until  some  items  fail, 
then  the  net  attrition  loss  just  described  must  be  increased  by  the  number 
of  units  required  to  "charge"  the  repair  pipeline  assuming  that  TAT  is 
shorter  than  the  PCLT  or  PCLT  plus  one  quarter.   The  respective  formulas 
for  demand  of  repairables  are  therefore: 


PCTT  TAT 

D„nT„   =   D.*[— ; — *(i  -  CRR*RSR)  +  -V--*CRR*RSRj 


D„-Tmtl   =   DA*[*PCyr+1'*(l  -  CRR*RSR)  +  ^*CRR*RSR] 

PCjLl-rl         A        4  4 
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APPENDIX  B:   THE  COSDIF  FORMULA 

B.l  THE  PHILOSOPHY 

The  COSDIF  formula  was  provided  by  Alan  Kaplan,  of  the  Army's  Inventory 
Research  Office,  who  was  a  member  of  the  team  that  developed  the  DODINST 
4140.42  [7].   It  is  derived  from  the  Army's  wholesale  range  model.   It  can  be 
viewed  as  the  results  of  a  simple  decision  model  as  illustrated  in  Figure  B.l 


Decision 


States  of  Nature 


No  demand  during  DDP 


Demand  during  DPP 


Make  a 

Costs  of  procurement 

Two  years  of  average 

provisioning 

plus  two  years '  holding 

annual  variable  costs 

buy 

costs 

Do  not  make  a 

Costs  of  spot  buys 

provisioning 

No  costs 

during  first  year  and 

buy 

average  annaul  variable 
costs  for  second  year 

Figure  B.l:   A  Provisioning  Decision  Matrix 

The  expected  costs  of  each  decision  can  be  evaluated  when  all  the  costs 
and  the  probability  values  for  the  states  of  nature  are  known.   That  decision 
which  corresponds  to  the  least  expected  cost  is  then  the  optimal  one. 

The  COSDIF  formula  given  below  determines  the  difference  between  the 
expected  costs  of  making  a  buy  and  not  making  a  buy. 

If  the  COSDIF  value  is  negative,  then  the  costs  associated  with  not  making 
a  buy  are  greater  than  those  for  making  the  buy,  and  is  therefore  optimal  to 
make  the  buy.   If  the  COSDIF  is  positive,  then  making  a  buy  is  more  costly. 
If  the  COSDIF  is  zero,  either  decision  would  be  optimal  but  not  making  a  buy 
is  less  work  for  the  provisicner  so  no  buy  is  made. 
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COSDIF  =   (D  Id)  *  {CP  +  2*H*C*(R+D__/FP)} 
o   T  SS 


+  (1-D  |D  )  *  {CP*FP  +  H*C*DSS/(2*FP)  +  D  *Cl) 


-  (1-dJd  )  *  {DT*(CSP  +  PLT*X/4)  +  DT*C*P> 


where: 


D  -   =   steady  state  annual  demand  (attrition) 

DT   =   total  steady  state  annual  demand  (TWAMP   *ERF) 

D  I D^   =   probability  of  no  demand  in  two  years,  given  an  annual 
steady  state  demand  forecast  (total)  of   D  .   This 
value  will  be  obtained  from  the  conditional  probabili- 
ties stored  in  th< 
rounded  up  at   .5, 


ties  stored  in  the  SCA  after  the  computed   D   is 


CP   =   cost  of  procurement  [the  small  purchase  value  is  used 
if  the  value  of  the  buy  quantity  is  less  than  or  equal 
to  the  small/large  purchase  breakpoint;  if  not,  the 
large  purchase  value  is  used.   The  specific  values  to 
be  compared  to  the  breakpoint  are  the  values  of  the 
buy  quantity   (R  +  D-q/FP)   and  the  economic  order 
quantity   (Dgg/FP)   for  lines  1  and  2,   respectively, 
of  the  formula] . 

H  =  holding  cost  rate 

C  =  unit  price 

R  =  reorder  point  quantity 

FP  =  frequency  of  procurement 

CI  =  cost  of  issuing  stock 

CSP  =  cost  of  spot  procurement 

:'..'  =  production  lead  time 


S 


implied  shortage  cost  [the  larger  of  the  shortage  cost 
for  CCSDIF  (specified  in  the  SCA  cog  constants)  and 
the  cost  to  hold,   H*C] 


P   =   spot  buy  premium  rate 

In  examining  the  COSPIF  formula  we  see  that  the  second  and  third  lines 
correspond  to  the  outcomes  of  the  alternative  decisions  when  chere  is  demand 
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during  the  DDP  (see  the  last  column  of  Figure  B.l).   However,  they  show  only 
the  average  annual  variable  costs  for  one  year  and  the  expected  costs  of 
making  spot  buys  for  one  year.   This  is  because  the  average  annual  variable 
costs  for  the  second  year  cancel  each  other  when  the  cost  difference  is 
taken  according  to  [7]. 

B.2   A  CRITIQUE 

The  COSDIF  formula  serves  as  the  range  model  in  the  development  of  the 
provisioning  budget  constraint.   As  a  consequence  it  is  appropriate  to  examine 
its  elements  and  to  critique  the  assumptions.   The  purpose  of  this  critique 
is  to  address  certain  details  of  the  formula  so  that  the  reader  will  be  able 
to  identify  where  improvements  in  the  current  use  of  the  model  by  the  Navy 
can  be  made.   It  is  not  possible  to  change  the  formula  without  approval,  and 
recent  attempts  to  do  so  have  not  been  successful.   These  will  be  discussed 
below. 

A  rather  puzzling  assumption  involves  the  use  of  the  steady-state  demand. 

It  is  especially  so  when  we  realize  that  the  depth  model  does  not  use  the 

steady  state  value.   Perhaps  it  is  true  that,  because  cose  differences  are 

being  used,  use  of   D    instead  of   D,   should  make  little  relative  difference, 

Sb  A 

This  should  be  evaluated. 

The  COSDIF  model  considers  only  two  states  of  nature  with  respect  to  the 
occurrence  of  demands  over  the  two  year  demand  development  period — zero  demands 
or  positive  demands.   SPCC  uses  probability  estimates  of  these  states  which 
are  a  function  of  the  expected  mean  demands.   The  probabilities  are  given  in 
Reference  [10] .   The  probabilities  used  by  SPCC  are  more  appropriate  than  those 
in  [7]  since  the  latter  were  based  on  inventory  records  for  non-Navy  items. 
In  fact,  [7]  required  each  service  to  come  up  with  better  values  within  two 
years.   The  probabilities  used  by  SPCC  can  be  further  improved  by  developing 
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values  for  each  two-digit  cog  as  was  attempted  in  [10] .   This  is  strongly 
recommended  since  cog  classes  do  represent  items  of  a  similar  nature.   It  is 
interesting  to  note  that  ASO  has  apparently  continued  to  use  the  probability 
values  given  in  [7]  for  all  of  its  cog  classes  [6]. 

It  could  be  argued  that  greater  resolution  in  the  demand  distribution 
should  be  considered  in  the  COSDIF  expression — perhaps  estimates  of  individual 
positive  demand  quantities  for,  say,  quarterly  or  yearly  periods  over  the 
two-year  DDP .   Even  if  the  probability  was  kept  at  a  single  value,  the  form 
of  the  average  annual  variable  costs  to  stock  an  item  over  the  two  years  sug- 
gests that  several  buys  will  be  made  and  orders  received.   The  buys  might  be 
made  but,  as  was  noted  in  the  main  part  of  this  report,  no  deliveries  are  ex- 
pected until  after  the  end  of  the  DDP.   In  fact,  the  average  variable  cost 
expressions  assume  that  steady  state  is  reached  instantaneously  at  POC  (as 
the  use  of   D  „   implied  earlier) .   A  switch  from  a  no-stock  to  a  stocking 
mode  between  the  first  and  second  year  is  virtually  impossible  because  of  the 
leadtimes.   The  spot  buy  situation  will  continue  through  all  of  the  DDP  under 
current  replenishment  leadtime  values.   There  will  be  no  holding  costs  since 
no  inventory  will  be  on  hand. 

The  term  of  the  COSDIF  equation  corresponding  to  the  zero-demand  state 
of  nature  is  the  closest  to  being  correct  ever  the  two  year  DDP.   The  impact 
o£   stocking  and  holding  for  two  years  is  significant  and  is  not  unexpected  in 
light  of  the  concern  in  [7]  about  excessive  unused  stocks  of  provisioned  items. 

Considering  the  spot  buy  term  in  more  detail,  it  is  interesting  to  note 
that  there  appears  to  be  no  explicit  term  for  an  item's  essentiality.   In  an 
earlier  document  (DODINST  4140.39  [6]),  the  replenishment  model  included  such 
a  term.   If  one  searches  through  [7]  for  a  reason  why  essentiality  was  not 
included,  the  only  consideration  that  is  found  is  in  Paragraph  IV. E.: 
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Non-demand-based  items,  which  do  not  meet  the  insurance  item  criteria 
for  wholesale  level  stockage  may  be  stocked  in  the  wholesale  system 
only  if  there  is  an  overriding  requirement  to  do  so  based  upon  their 
essentiality  in  a  selected  weapon  system.   In  this  instance,  an  item 
will  be  stocked  as  a  Numeric  Stockage  Objective  (NSO)  item. 

Reference  [7]  also  notes  that 

no  deviation  to  the  basic  formula  is  authorized,  the  use  of  cost  data 
more  relevant  to  a  particular  inventory  should  provide  sufficient 
flexibility  to  account  for  the  variability  required  in  support  of 
different  systems.   Data  developed  for  ICP  wholesale  level  systems 
in  accordance  with  DODINST  4140.39  should  be  useful  in  establishing  the 
most  appropriate  holding  costs,  procurement  costs,  and  implied  shortage 
cost . 

Additionally,  a  constraint  is  imposed  such  that  the  cost  of  a  requisition  on 

backorder  must  be  at  least  equal  to  the  cost  to  hold  an  equivalent  amount  of 

material.   It  may  also  be  desirable  to  use  a  larger  value  for  the  holding  cost 

in  the  COSDIF  than  is  used  in  Reference  [6]  to  reflect  a  higher  obsolescence 

probability  for  marginally  stocked  items  [7]. 

A  variation  in  the  COSDIF  which  includes  an  essentiality  term  was  presented 
in  Reference  [9].   However,  the  goal  of  that  reference  was  to  justify  the 
variable  threshold  model  (which  does  include  an  essentiality  term  in  the  RISK 
equation  used  to  determine  the  depth)  as  an  alternative  model  for  implementation 
of  the  budget  constraint.   It  is  important  to  note  that  D55,  which  is  SPCC's 
implementation  of  the  COSDIF,  does  not  use  this  variation  in  creating  the 
budget  constraint. 

In  July  1981,  the  Office  of  the  Chief  of  Naval  Operations  [4]  proposed 

changing  Reference  [~]  to  correspond  to  the  form  of  the  COSDIF  in  [9j.   The 

response  from  the  Deputy  Assistant  Secretary  of  Defense  (Logistics  and 

Material  Management)  was  [11] : 

While  we  understand  the  thrust  of  your  initiative  to  revise  DoD 
Instruction  4140.42,  we  do  not  agree  with  the  necessity  of  making 
the  requested  changes.   This  instruction  currently  recognizes  essen- 
tiality in  the  initial  spares  computation  and  allows  stockage  of  essen- 
tial items  even  if  they  do  not  qualify  for  stockage  through  the  use 
of  the  cost  differential  tables. 
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When  a  shortage  cost  is  considered  in  any  inventory  model  supporting  a 
weapon  system,  it  is  implicitly  a  measure  of  the  essentiality  of  the  item. 
The  value  that  can  be  used  for   A    in  the  COSDIF  equation  is  not  constrained 
by  References  [7]  and  [17]  to  be  that  of  the  replenishment  model.   In  fact, 
the  values  used  in  the  replenishment  models  are  tied  to  a  replenishment  budget 
constraint  in  any  given  fiscal  year.   As  a  consequence,  the  shortage  cost 
values  in  the  replenishment  model  can  change  each  year.   It  should  be  obvious 
then  that  such  values  are  inappropriate  for  provisioning. 

Control  over   X   will  have  a  significant  effect  on  the  development  of 
the  budget  constraint.   The  larger  the  A   value,  the  larger  will  be  the 
range  of  demand-based  items  used  in  developing  the  budget  constraint.   The 
bottom  line  is  to  obtain  a  larger  budget  constraint.   The  determination  of 
better  A   values  should  be  viewed  as  an  important  first  step  of  improving 
the  provisioning  process.   The  ASO  provisioning  model  attempts  to  do  this. 
It  adjusts   A    in  an  effort  to  gain  an  SMA  of  85%  for  the  depth  allowed  by 
[7]  for  determining  the  budget  constraint.   This  is  done  separately  for   1R 
cog  and  for  the  group  of   2R  and   6R  cog  items.   The  ASO  model  develops  a 
cost  constraint  for  five  different  values  of   A    using  application/operation 
D54  [17].   For  each  A  ,   the  depths  of  the  items  are  used  to  estimate  the 
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SMA  which  would  result.   A  curve  is  then  constructed  of  SMA  versus  investment. 

Additional  A   values  as  needed  are  used  to  generate  additional  points  as  the 
s 

process  zeros  in  en  the  85%  SMA  point.   The  investment  value  for  the  85%  SMA 
point  is  then  established  as  the  budget  constraint.   The  actual  range  and 
depth  of  provisioned  items  are  computed  using  D52  [18] .   D52  uses  a  load  list 
approach  developed  prior  to  Reference  [7];  the  goal  of  that  approach  is  to 
minimize  expected  units  short.   The  A   value  generated  bv  D54  is  used  only 
tO  generate  the  budget  constraint.   D52  uses  a  Lagrange  multiplier  approach  to 
implement  thac  budget. 
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It  is  debatable  whether  SMA  is  the  appropriate  effectiveness  goal  to  use 
in  light  of  the  current  trend  toward  A  .   Additionally,  SMA  is  measured  only 
for  those  items  passing  the  COSDIF  range  criterion.   However,  the  approach 
used  in  the  ASO  model  is  certainly  worth  considering  with  availability  goals 
used  in  place  of  the  SMA  goal. 
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APPENDIX  C 

This  appendix  gives  a  simple  counter-example  to  show  that  maximization 
of  "pseudo-availability"  as  defined  by  the  relationship 


n         MTBF. 

a     =    n 


o      .  MTBF.  +  MTTR.  -I-  MSRT.(s.) 
11       i       11 


does  not  yield  the  same  solution  to  the  stockage  allocation  problem  as 
minimization  of  mean  supply  response  time  (MSRT)  or  time-weighted  units  short 
(TWUS) .   Thus,  the  claim  that  "sparing  to  availability"  is  equivalent  to 
sparing  to  minimize  logistics  delay  time  is  not  true. 

Example:   Consider  a  system  consisting  of  two  components  having  the  parameter 
values  given  in  Table  CI: 

PARAMETER  1  2 


MTBF  0.200         0.100 

PCLT  1.0  1.0 

MTTR  0.0822        0.0274 

C  5.  10. 

TABLE  CI:   SYSTEM  PARAMETERS 

where  PCLT  is  the  procurement  leadtime  in  years;   C   is  the  unit  cost; 
MTBF  =  1/a   is  the  mean  time  between  failures;  and  MTTR  is  the  mean  time 
to  repair.   Suppose  ve  are  allowed  a  total  provisioning  budget  of   B  =  S20. 
Then  it  is  clear  chat  there  are  only  3  undominated  feasible  solutions  to  the 
allocation  problem:   (4,0),  (2,1),   and   (0,2).   Table  C2  gives  the  mean 
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supply  response  times  for  each  of  the  feasible  values  of  stockage  levels  for 
components  one  and  two. 

MSRT(s) 


s 

Item 

C 

) 

1 

> 

3 

4 

1 

182. 

50 

124 

.00 

79. 

51 

47.80 

26.83 

2 

182 

50 

149 

.65 

120 

45 

X 

X 

TABLE  C2:   MSRT  vs.  s 

From  the  information  in  Tables  CI  and  C2  it  is  easy  to  determine  the  overall 
mean  supply  response  time  and  the  system  availability  for  each  of  the  three 
candidate  solutions.   These  values  are  given  in  Table  C3. 


CANDIDATE  SOLUTION 


MSRT 


(4,0) 
(2,1) 
(0,2) 


130.61  0.0896  = 

126.27*  0.0744 

141.13  0.0559 


TABLE   C3:       SYSTEM  MSRT   AND  A 

o 

The  values  in  Table  C3  reveal  that  the  solution  (2,1)  is  optimal  with 
respect  to  total  MSRT  whereas  the  solution  (4,0)  is  optimal  with  respect 
to  system  availability! 
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